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ABSTRACT 

This paper studies the mean angular momentum associated with the collision of two 
celestial objects in the earliest stages of planet formation. Of primary concern is the 
scenario of two rarefied preplanetesimals (RPPs) in circular heliocentric orbits. The 
theoretical results are used to develop models of binary or multiple system formation 
from RPPs, and explain the observation that a greater fraction of binaries originated 
farther from the Sun. At the stage of RPPs, small-body satellites can form in two 
ways: a merger between RPPs can have two centers of contraction, or the formation of 
satellites from a disc around the primary or the secondary. Formation of the disc can 
be caused by that the angular momentum of the RPP formed by the merger is greater 
than the critical angular momentum for a solid body. One or several satellites of the 
primary (moving mainly in low-eccentricity orbits) can be formed from this disc at any 
separation less than the Hill radius. The first scenario can explain a system such as 2001 
QW322 where the two components have similar masses but are separated by a great 
distance. In general, any values for the eccentricity and inclination of the mutual orbit 
are possible. Among discovered binaries, the observed angular momenta are smaller 
than the typical angular momenta expected for identical RPPs having the same total 
mass as the discovered binary and encountering each other in circular heliocentric 
orbits. This suggests that the population of RPPs underwent some contraction before 
mergers became common. 

Key words: Solar system: formation; minor planets, asteroids; Kuiper belt; planets 
and satellites: formation 



1 INTRODUCTION 

1.1 The case for rarefied preplanetesimals 

From the 1950's to the 1980's, many authors (e.g. Safronov 
1969; Goldreich & Ward 1973) believed that planetesimals 
originated from rarefied dust condensations (RDCs): over- 
dense regions of the disc that are gravitationally bound, and 
may contain solid particles or bodies, but are still gaseous 
or diffuse. During this period, the mechanisms behind RDC 
formation and the coagulation of larger secondary conden- 
sations were studied by several authors (e.g. Safronov 1969, 
1995; Eneev & Kozlov 1981, 1982; Pechernikova & Vityazev 
1988). The duration of the condensation stage was thought 
to be longer at greater distances from the Sun, where the 
RDCs have lower densities. According to an early model by 
Safronov (1967, 1969), the time required for RDCs to form 
into solid bodies was about 10 4 yr in the terrestrial zone 
and 10 6 yr at the distance of Jupiter. The corresponding 
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increases in mass were 10 2 and 10 3 from the initial values. 
Safronov concluded that while the rotation of initial RDCs 
impedes contraction, they also become denser when they 
combine. 

In the models of RDC contraction considered by Myas- 
nikov & Titarenko (1989, 1990), collisions between conden- 
sations were not taken into account. The times needed for 
RDCs to become solid planetesimals could exceed several 
million years, depending on the optical properties of the 
dust and gas and the type and concentration of short-lived 
radioactive isotopes in RDCs. Depending on the concentra- 
tion, the times varied by a factor of more than 10 and could 
exceed 10 Myr at the concentration greater than 0.02. 

In the 1990's, scenarios of planetesimal formation 
through the gravitational instability of RDCs were frus- 
trated by the phenomenon of self-generated mid-plane tur- 
bulence. Models involving the hierarchical accretion of plan- 
etesimals from smaller solid objects became more popular 
(e.g. Weidenschilling 2003). In recent years, however, new 
arguments (e.g. Makalkin & Ziglina 2004; Johansen et al. 
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2007; Cuzzi, Hogan & Shariff 2008; Lyra et al. 2008; Jo- 
hansen, Youdin & Mac Low 2009) have arisen in favour of 
RDCs, or as they can be also called rarefied preplanetesi- 
mals (RPPs) (condensations, clusters, and clumps). Unlike 
earlier models, however, the overdense regions could include 
metre-sized boulders and smaller solid bodies. Some results 
from these more recent models are summarised below. 

Makalkin & Ziglina (2004) showed that after the sub- 
disc reaches a certain critical density, its inner, equatorial 
layer could be gravitationally unstable. This is possible be- 
cause the inner disc, unlike the two subsurface layers, con- 
tains no shear turbulence. In their model, RDCs form trans- 
Neptunian objects (TNOs) with diameters up to 1000 km 
in ~10 6 yr. 

Johansen et al. (2007) found that gravitationally bound 
clusters can form with masses comparable to dwarf plan- 
ets. They determined that peak densities in the clusters ap- 
proached W 4 p g as (where p gas is the average gas density) 
after only seven orbits. Lyra et al. (2008) noted that strong 
drag forces in the disc might delay gravitational collapse 
of the clusters. Cuzzi et al. (2008) showed how self-gravity 
could stabilise dense clumps of millimetre-sized particles. 
Such clumps form naturally in 3D turbulence, forming dif- 
fuse but cohesive sandpiles on the order of 10-100 km in di- 
ameter. For chondrules of radius r c — 300/im, they obtained 
a 'sedimentation' timescale of about 30-300 orbital periods 
at 2.5 au. (The timescale is proportional to rj 1 .) Youdin 
(2008) summarised various mechanisms for particle concen- 
tration in gas discs, including turbulent pressure maxima, 
drag instabilities, and long-lived anticylonic vortices. 

Ida, Guillot & Morbidelli (2008) studied the accretion 
and destruction of solid planetesimals in turbulent discs, and 
found that accretion proceeds only for planetesimals with di- 
ameters above 300 km at 1 au and above 1000 km at 5 au. 
Their analytical arguments were based on fluid dynamical 
simulations and orbital integrations. They concluded that 
some mechanism must be capable of producing Ceres-mass 
planetesimals on very short timescales. Based on the ob- 
served size-frequency distribution and collisional evolution 
of asteroids, Morbidelli et al. (2009) concluded that the ini- 
tial planetesimals had to range from one hundred to several 
hundreds of km in size, probably up to 1,000 km. 

Prior to these new arguments in favour of RPPs, I (Ipa- 
tov 2001, 2004) had suggested that TNOs with mildly ec- 
centric orbits (e < 0.3) and diameters greater than 100 km 
(such as Pluto and Charon) could result from the compres- 
sion of large RPPs with semimajor axes a > 30 au rather 
than the accretion of small, solid planetesimals. I also pro- 
posed that some planetesimals with diameters d ~ 100-1000 
km in the feeding zone of the giant planets, some planetes- 
imals with d ~ 100 km in the terrestrial zone, and some 
large main-belt asteroids could have formed by direct com- 
pression. Some smaller objects (TNOs, planetesimals, and 
asteroids) with d < 100 km could be debris from large ob- 
jects, while others could have formed by direct compression 
of preplanetesimals. 

1.2 Formation of binaries 

The frequency of binary systems and the ratio of sec- 
ondary mass to primary mass are greater among classi- 
cal TNOs with inclinations i < 5° than among main- 



belt asteroids (Richardson & Walsh 2006; Knoll 2006; 
www.johnstonsarchive.net / astro / asteroidmoons.html) . 

Planetesimal formation models based on the accretion 
of solid bodies have put forth several hypotheses on the for- 
mation of binaries. For example, several papers (e.g. Dores- 
soundiram et al. 1997; Durda et al. 2004; Canup 2005; Walsh 
2009) have been devoted to the mechanism of catastrophic 
collisions. Weidenschilling (2002) studied the collision of two 
planetesimals within the sphere of influence of a third body. 
Goldreich, Lithwick & Sari (2002) proposed that the pri- 
mary body could capture a secondary component passing 
inside its Hill sphere due to dynamical friction from sur- 
rounding small bodies or gravitational scattering by a third 
large body. Funato et al. (2004) studied a model in which 
a low-mass secondary component is ejected and replaced by 
a third body in a wide but eccentric orbit. Astakhov, Lee 
& Farrelly (2005) studied binary system formation in four- 
body simulations with solar tidal effects. Gorkavyi (2008) 
proposed a multi-impact model. Six mechanisms of binary 
formation were discussed by Doressoundiram (1997). More 
references can be found in papers by Richardson & Walsh 

(2006) , Noll (2006), Petit et al. (2008), Noll et al. (2008), 
and Scheeres (2009). 

We agree with Cuk (2007), Pravec, Harris & Warner 

(2007) , Walsh, Richardson & Michel (2008), and several 
other authors that binary systems with small primaries 
(such as near-Earth objects) arise mainly from the rota- 
tional breakup of 'rubble piles', for example due to increas- 
ing spin from the Yarkovsky-O'Keefe-Radzievskii-Paddack 
(YORP) effect. However, I believe that collisions could also 
push small bodies beyond their critical spin limits. 

Some scientists believe that the formation of small-body 
binaries could be similar to the formation of binary stars 
within fragmented discs (Alexander, Armitage & Cuadra 
2008). Nesvorny (2008) noted that an excess of angular mo- 
mentum prevents the agglomeration of all available mass 
into solitary objects during gravitational collapse. 

1.3 Problems studied in this paper 

In contrast with models of binary formation that consider 
solid bodies, in 2004 I proposed that a considerable frac- 
tion of trans-Neptunian binaries (including Pluto-Charon) 
could form while rarefied preplanetesimals moving in almost 
circular orbits undergo compression (Ipatov 2004). In 2009, 
I set forth two detailed models of binary formation during 
this stage (Ipatov 2009). They are repeated here and further 
discussed in Sections 2.2 and 2.3. These models can explain 
trans-Neptunian binaries, and some asteroid binaries with 
large primaries (with diameters of at least 50 km). Small 
asteroid binaries could mainly arise from collisions between 
larger solid bodies. 

The problem of binary formation from RPPs has sev- 
eral aspects (e.g. simulations of contraction, collisions be- 
tween preplanetesimals). In the present paper, I pay par- 
ticular attention to calculating the angular momentum of 
two colliding RPPs about their centre of mass (Section 3) . I 
also compare the predicted angular momenta to the observed 
momenta of discovered binary systems (Section 4). Other as- 
pects of binary formation at the preplanetesimal stage (e.g. 
mergers of colliding preplanetesimals) are discussed briefly 
in Section 5, and may be the subject of future publications. 
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We show that if RPPs existed during the period permit- 
ting encounters to occur between preplanetesimals within 
their Hill spheres, then it is possible to explain the differ- 
ences between TNO and asteroid binaries as well as some 
of the peculiarities of observed binaries. The material pre- 
sented in Section 3 also helps us better understand the angu- 
lar momenta of planetesimals, trans-Neptunian objects, and 
asteroids, although these momenta could change consider- 
ably through gravitational interactions and collisions after 
all bodies have solidified. 



2 SCENARIOS OF BINARY FORMATION 
DURING THE RAREFIED 
PREPLANETESIMAL STAGE 

2.1 Application of previous solid-body scenarios 
to preplanetesimals 

At present, most small bodies move in orbits with eccen- 
tricities exceeding 0.05. Solid planetesimals would have be- 
gun with much smaller eccentricities, but mutual gravitation 
tends to increase the eccentricities and inclinations of their 
orbits (e.g. fpatov 1988, 2007). RPPs, on the other hand, 
would remain in almost circular orbits with very small in- 
clinations before the moment of formation of a binary. The 
smaller growth of RPP orbital eccentricities could be be- 
cause the stage of preplanetesimals was short, and circula- 
tion of an orbit of an object (e.g. a preplanetesimal or a 
planetesimal) of a fixed mass due to the influence of gas and 
dust could be greater for a larger size object. Due to their 
larger sizes, preplanetesimals collided more often than solid 
planetesimals per unit of time. 

If we consider almost circular heliocentric orbits, the 
typical minimum distance between the centres of mass of 
the objects could be smaller than that between two objects 
entering the Hill sphere from eccentric heliocentric orbits. 
Heliocentric orbits with high eccentricities more often give 
rise to hyperbolic relative motion. In contrast, for nearly 
circular orbits, the trajectory of relative motion inside the 
sphere can be complicated (e.g. Ipatov 1987; Greenzweig & 
Lissauer 1990; fwasaki & Ohtsuki 2007). In this case, the 
objects can move inside the sphere for a long time, and the 
relative distance can change considerably. The role of tidal 
forces could also be greater for rarefied objects. 

The above discussion supports my belief that models 
of binary formation due to gravitational interaction or colli- 
sions occurring inside the Hill sphere of the future primary, 
which have been considered by several authors for solid ob- 
jects, could be even more effective for rarefied preplanetesi- 
mals. 

2.2 Binary formation from two centres of 
contraction 

Rarified preplanetesimals contract over time, and may ex- 
hibit centrally concentrated radial density functions. It is 
reasonable to suppose that in some cases, the collision of 
two such objects results in a rotating system with two cen- 
tres of contraction. The end result of gravitational collapse 
would be a binary system with almost the same total mass 
and components separated by a large distance (such as 2001 



QW322). For such a scenario, the values of the eccentricity of 
the present mutual orbit of components can have any values 
less than 1. 

If the original RPPs were similar in size to their Hill 
spheres, the separation L of the solid binary resulting from 
the combined system could range up to the radius m of the 
Hill sphere. More often, however, the preplanetesimals were 
much smaller than their Hill spheres. The separation dis- 
tance L would then be less than m ■ For discovered binaries 
with L/r p > 100, the ratio r a /r p of the radii of secondary 
to primary components is greater than 0.5. This fact may 
imply that when the masses of two rarified preplanetesimals 
differ greatly, only one centre of contraction survives. Note 
that other authors (e.g. Descamps & Marchis 2008) have 
suggested different explanations for the above ratio. 

2.3 Binary formation due to excessive angular 
momentum 

In my opinion, some binaries could have formed from RPPs 
that obtained more angular momentum than the maximum 
possible value for a solid body. Let v a denote the velocity 
of a particle on the surface of a rotating object, and let 
v c f be the minimum velocity at which a particle can leave 
the surface. As a rotating rarefied preplanetesimal contracts, 
some of the material with v a > v c f could form a cloud 
(that transforms into a disc) that moves around the pri- 
mary. For a spherical object, v a is greatest at the equator 
and v c f can be close to the circular velocity. For relatively 
condense preplanetesimals, their collision can eject material 
into a satellite- forming disc. The disc can capture dust and 
boulders that enter the Hill sphere after the encounter of 
RPPs. 

One or several satellites of the primary could form in 
this cloud (similar to typical models of formation of satellites 
of planets; see Woolson 2004, for example). The eccentrici- 
ties of such satellites would usually be small. As we show at 
the end of Section 3.3, this scenario allows the formation of 
satellites at practically any distance from the primary up to 
the Hill radius. The initial radius of the 'parent' preplanetes- 
imal, however, could be greater than the separation between 
components of the formed binary. 

2.4 Hybrid scenario and the formation of 
elongated bodies 

The two scenarios presented in Sections 2.2 and 2.3 could 
take place at the same time. In addition to massive primary 
and secondary components, smaller satellites moving around 
either body could be formed. For binaries formed in such a 
way, the separation distance between the main components 
can be large or small. 

It is possible that massive yet rarified binary compo- 
nents merge when their central parts became dense enough. 
In such a case, the form of the solid body obtained could 
differ greatly from an ellipsoid. For example, this situation 
could give rise to an elongated, bone-like body where both of 
the original components are visible. As in other scenarios, 
the solid body may have one or more small satellites. For 
example, (216) Kleopatra (with dimensions of 217 x 94 x 81 
km) could have formed in this manner. 
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Even more rarely, several partly compressed compo- 
nents could merge simultaneously. In this case the form of 
the solid body would be complicated, to some extent 're- 
membering' the forms of its components. 



3 THE ANGULAR MOMENTUM OF 
COLLIDED RAREFIED 
PREPLANETESIMALS 

3.1 Basic model of an encounter between 
preplanetesimals 

In order to analyse the scenarios discussed in Section 2, Sec- 
tions 3.2-3.4 will study the angular momentum K s of two 
colliding RPPs relative to their centre of mass and the pe- 
riod T s of the solid primary formed by contraction. Maxi- 
mum and typical values of the separation distance L between 
binary components for the above scenarios are discussed in 
Section 3.3. 

Most of the formulae and calculations below are pre- 
sented for a 'basic' model: the merger of two spherical RPPs 
that before the collision were moving in coplanar, circular, 
heliocentric orbits. The mass of the resulting object is equal 
to the sum of the two preplanetesimal masses. The rotating, 
merged system then contracts and tranforms into a solid 
planetesimal. The angular momenta presented in Table 1 
were obtained assuming that the radii of the two preplan- 
etesimals are equal to their Hill radii. Based on the 'basic' 
model, we discuss more complicated models. 

Scientists who study the formation of preplanetesimals 
generally consider almost circular orbits and very small incli- 
nations. These conditions minimise the relative velocities of 
neighbouring preplanetesimals. The probability of a merger 
is smaller for greater eccentricities. In all existing models of 
planetesimal discs, the mean inclination (in radians) of the 
orbits is assumed to be smaller than the mean eccentricity. 
Relaxing this condition on the inclinations will not change 
the order of magnitude of the angular momentum of two 
colliding preplanetesimals, but can change its direction (see 
Sections 4.3 and 5.3). 

For formation of a binary, it is not necessary that every 
encounter between preplanetesimals (i.e., when one object's 
centre of mass passes through the Hill sphere of the other) 
results in a merger. It is enough that for some preplanetesi- 
mals, at least one such encounter occurred during their life- 
times. For almost circular heliocentric orbits, the minimum 
distance between encountered material points with masses of 
preplanetesimals can be several orders of magnitude smaller 
than the Hill radius rn, and some of the material points 
can move inside the sphere for a long time. For example, in 
an earlier work (Ipatov 1987) I presented a case where two 
planetesimals move together for more than half an orbital 
period (see also discussions in Sections 2.1 and 5.1). There- 
fore, the values of K s and T s obtained for the 'basic' model 
can also be true for encounters between preplanetesimals 
much smaller than their Hill spheres. The sizes and mergers 
of encountered preplanetesimals are discussed in Section 5.1. 

Not all boulders could be captured during an encounter, 
and some material could be lost at the stage of contraction. 
Therefore, the mass of the formed binary (or a preplanetes- 
imal) may be smaller than the sum of the masses of the 



two colliding preplanetesimals. Likewise, the final value of 
K s may be smaller than the total angular momentum of the 
RPPs. 

Sections 3.2 and 3.3 present formulae for T s and K s 
after a single collision in the basic model. For multiple col- 
lisions, the angular momenta are summed. In my opinion 
(see Section 4.3), most of the angular momentum possessed 
by a trans-Neptunian binary usually resulted from a single 
collision. 

Naturally, a present-day binary can differ significantly 
from the primordial binary formed by contracting RPPs. For 
example, tidal forces can increase the separation between 
binary components (e.g. the Pluto-Charon and Antiope- 
S/2000(90) systems) and decrease the period of axial rota- 
tion (as occurred for the Earth-Moon system) on time scales 
shorter than the age of the Solar system. Although a given 
primordial binary system might have been more compact 
than its corresponding discovered system, the total angular 
momentum could remain the same. The angular momenta 
of individual components could also change due to collisions 
with other small bodies and the YORP effect. Finally, the 
angular momentum of a contracting preplanetesimal could 
decrease due to tidal interactions with the Sun. 

It is likely that most primordial planetesimals were el- 
lipsoidal. The velocities of particles close to the equator of a 
rotating contracting preplanetesimal were greater than those 
of particles close to the poles, and more particles could leave 
the equator. Therefore, some contracting preplanetesimals 
could be elongated along the axis of rotation. On the other 
hand, a surrounding disc could feed material to the equa- 
tor, resulting in an oblate shape. In either case, the equa- 
tor would be approximately circular. In the 'rubble pile' 
model considered by Walsh et al. (2008), the axial radius 
was smaller than the equatorial radius. The formation of 
highly elongated primordial small bodies was discussed in 
Section 2.4. 

3.2 Principal formulae for calculation of the 
angular momentum of two colliding 
preplanetesimals 

Previous papers (e.g. Lissauer et al. 1997; Ohtsuki & Ida 
1998) devoted to the formation of axial rotation of celestial 
bodies have mainly studied the final rotation rates of planets 
accreted from a disc of solid planetesimals. My own studies 
on the formation of axial rotation of rarefied preplanetesi- 
mals, solid planetesimals and planets were presented in two 
preprints (N 101 and N 102) published in Russian in 1980 by 
the Institute of Applied Mathematics of the USSR Academy 
of Sciences. These studies are summarised in Ipatov (2000) , 
which cites many papers devoted to the formation of spin. 
The preprints focussed mainly on the axial rotation of plan- 
ets, but the formulae used are applicable to a wide variety of 
cases. Sections 3.2-3.4 present similar formulae. In Sections 
4-5, the formulae are applied to the studies of the formation 
of binaries. 

Consider an inelastic collision between two non-rotating 
spherical objects (preplanetesimals or planetesimals) with 
masses mi and rri2 and radii ri and T2, moving around the 
Sun in coplanar circular orbits. The first object is closer to 
the Sun. Its semimajor axis is denoted by a; its circular ve- 
locity can be calculated from the relationship vi = v c — 
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(G Ms/a) 1 / 2 , where G is the gravitational constant and Ms 
is the mass of the Sun. The angular momentum of the col- 
liding system relative to its centre of mass is 



v T (ri + r2)mim2/(mi + 7772), 



(1) 



where v T is the tangential component of the velocity (with 
respect to the centre of mass of the first object) at which the 
objects collide. If the difference Aa in the orbits' semimajor 
axes equals 0(r x + r 2 ), we have 

v T = kev c (n +r 2 )/a = ke(G Ms) 1/2 {ri + r 2 )a~ 3/2 
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(2) 



(3) 



(4) 



and 

K s ee K. e = ke(GM s ) 1/2 ( ri +r 2 ) : . 

{mi + m2)a i i 2 - 

For r a — (n + r 2 )/a -C 0, one can obtain 

fc e w (1 - 1.50 2 ) 

and 

v rel « (l-O.750 2 ) 1/2 v c r a , 

where v re i is the absolute value of the relative velocity be- 
tween the two colliding objects. 

Below I derive equations (3) and (4). The angle between 
the lines connecting the Sun to the two encountering objects 
is a ~ sina ~ r a (l — 2 ) 1 ' 2 . Consider the triangle ABC, 
where points A and B are the positions of the objects and 
point C is the intersection between the orbit of the first 
(faster) object and the line from the Sun to the second object 
(point B). The angle CAB is close to the angle 7 between 
the velocity of the first object and the direction from the 
first object to the second object (located farther from the 
Sun). We therefore consider sin 7 « 6. Taking into account 
that v c oc aT 1 ^ 2 and (1 + r a 6) _1 ^ 2 « 1 — r a Q/2, we obtain 
that the velocity of the second object equals v 2 ~ «i(l — 
r a 0/2). For the considered model, cos a « 1 and cos 7 « 
(1 — 9 2 ) 1//2 . We therefore obtain sin(7 + o?) « sin7 + acos7, 
and the difference between the tangential velocities of the 
two encountering objects is 

v T i — v T 2 = vi sin 7 — V2 sin(7 + a) 

« tji6 - wi(l - r a ©/2)(0 + (r„(l - © 2 ) 1/2 )(1 - 2 ) 1/2 ) 
«vir„(1.50 2 -l). 

This difference is valid for clockwise motions (negative an- 
gular velocities). Therefore, k& « (1 — 1.59 2 ). 

Let us take into account that the angle between ve- 
locities vi and 7J2 is equal to a, v 2 ~ vi(l — r a ©/2), 
sina w r a (l — O 2 ) 1 ^ 2 , and cosa « 1. Neglecting the term 
with r„, we obtain 

^rei ~ (vi — ^2 cos a) 2 + (v 2 sin a) 2 

« (9 2 /4 + (1 - 2 ))(« c r a ) 2 = (1 - 0.756 2 )(7; c r a ) 2 . 

Though the above derivation of equation (3) is not valid 
for small O, Below we can use it to estimate the mean values 
of fee, K s , and v T , by evaluating the integrals over 0. The 
reviewer noted that equations (3) and (4) are actually ac- 
curate even for = 0, as planetesimals 'colliding' (or more 
accurately, touching) on the same circular heliocentric or- 
bit have nonzero angular momentum relative to each other 
in the inertial frame despite being stationary in the rotat- 
ing frame. Only at larger does the Keplerian shear begin 



to dominate, resulting in a negative (i.e. retrograde) total 
angular momentum. 

The overall mean value of | fee | is 0.6. Among posi- 
tive values of fee the mean is |, and among negative values 
it is —0.24. The values of K s q and fee are positive in the 
range < < | X ^ 2 w 0.8165 and negative in the range 
0.8165 < < 1. The same ranges of for positive and neg- 
ative values of K s q were used by Eneev Kozlov (1982), who 
obtained them from a numerical integration of the three- 
body problem. The minimum value of fee is —0.5, and its 
maximum value is 1. If is uniformly distributed for ob- 
jects moving in the same plane, the probability of a single 
collision producing negative rotation is about i. The ratio 
of the sum of positive values of K s to the sum of all | K s | 
equals 0.925. This ratio, obtained by evaluating the integrals 
given above, is close to the ratio obtained in numerical ex- 
periments with initially circular orbits. The results of the 
experiments were presented in the aforementioned preprints 
in 1981. 

In the basic model, a new object forms in any encounter, 
and heliocentric orbits are considered to be circular at the 
time of encounter up to rjj. This model could work for some 
encounters of very rarefied preplanetesimals (see discussion 
in Section 5.1). Compared to the above model, the fraction of 
collisions with positive angular momentum would be smaller 
if one considers collisions of solid (or at least higher-density) 
objects and/or supposes that their orbits are circular only so 
long as the objects are widely separated (that is, the orbits 
deviate from a heliocentric circle when separation equals to 
rsr). Under the right conditions, it is even possible for the 
merged planetesimal to have negative momentum. There- 
fore, the dominance of positive angular momentum in the 
basic model does not contradict the negative spins obtained 
by authors who have used other models and considered col- 
lisions of solid bodies (e.g. Giuli 1968). 

Below we assume that a new spherical object of radius 
r/ and mass m/ = 7711+7712 forms as the result of a collision, 
and then contracts to radius 77 = r/ /fe r . If the angular 
momentum K s q relative to the centre of mass is the same 
as that of the system at the time of collision, then equation 
(2) yields the period of axial rotation: 



_ 2tt J s , (mi + m 2 ) 2 r 2 f a 3/2 

Js = —7^ = KT- 



K s 



77117772 (t*i + 72) 2 



(5) 



where kr = 0.8n xkQ 1 (G Ms)' 1 ^ 2 ■ The moment of inertia 
is J s = OAxmfrp where \ — 1 f° r a homogeneous sphere. 
At rj a — rf + r\, mi = 77i 2 , and 7*1 = r 2 , using equations 
(1) and (5), we then obtain 

T s = 2 2/3 fe T a 3/2 fc T : 2 « 1.67r X r 2 /(7j T ri) w 6.33xr//(i; T fc r ).(6) 
For mi ^> 7712 and 77 « ri/k r , we have 



kTira\ I m 2 )d i ^ 2 k r 2 w 0.87T xr\mi/(m,2V T k 2 ). 



(7) 



In the case of heliocentric orbits with eccentricity e, we can 
use v T « 2~ 1//2 v c e. 

The angular momentum of binary components (with 
masses m p and m s and radii r p and r s ) relative to their 
centre of mass equals 

Kscm = v b Lm p ms/(m p + m s ), (8) 
where L is the mean distance between the components (i.e. 
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the semimajor axis of their mutual orbit) and vt = (G(m p + 
ras)///) 1 / 2 is the mean velocity of their relative motion. The 
ratio of the momentum K a e obtained at the collision of two 
spheres - preplanetesimas to K scm is 



k K = 



K s@ _ k (M s ) 1/2 (m p + m s ) 1/2 m 1 m 2 {r 1 + r 2 f 



m p m s (mi + m^a?/ 2 !^ 1 / 2 



.(9) 



3.3 The spins of colliding preplanetesimals that 
have radii proportional to their Hill radii 

This subsection calculates K s q under the assumption that 
the radius of each preplanetesimal is proportional to the ra- 
dius of its Hill sphere, i.e. n = k H a(m i /'iMs) llz (i = 1, 2). 
The coefficient ku is arbitrary; as discussed in Section 3.1, 
we often consider ku = 1. The radii used to calculate K s q 
can be much larger than the physical radii of the preplan- 
etesimals. 

Let us consider mi = m 2 = m p — m s = m/ /2 = |-7rr 3 p 
and r p = r a = r, where p is the density of solid binary 
components of radius r. In this case, formulae (2), (8), and 
(9) yield 

K s@ = e-WklkeMg^G 1 ^ 1 ' 2 ™*/ 3 



(10) 



and 

k 2 K — 6Msk%kfr /it p L a 3 



-1/3 



arL- 1 p Ll ' i k%k\ i . (11) 



From equation (11) we obtain 
L ?s 1.5k%kjjk K 2 r H , 

where m ~ a(m//3Ms)'' 3 . For kK — ku = 1, we have L « 
0.5th at fee = 0.6, and L « r^r at fee = 0.816. Therefore, for 
an encounter between rarefied preplanetesimals moving in 
almost circular heliocentric orbits, a binary with components 
of equal mass could form with any separation up to the Hill 
radius. 

In Sections 4.2-4.3, 1 will compare the values of K s q and 
K scm for several discovered binaries. Using the dependence 
(11) of kK on p and L, one can estimate the uncertainty on 
kK for the binaries where these values are not well known. 

At r/ = kH P rTH, the ratio k r — rf /rf is proportional 
to a p 1 / 3 and does not depend on mass. It equals 133 at a = 1 
au, kn P r = 1, and p = 1 g cm~ 3 . For a solid planetesimal 
formed by contraction of a preplanetesimal of radius r/ , 
one obtains 

rn 3/2 , -2 -1/2 -2/3 

T s oc a ' k r oc a ' p ' . 

At x = lj fce = 0.6, a—1 au, fcn = kij P r = 1, and 
mi = m,2, the rotation period of a RPP formed from the 
collision of two Hill sphere-sized preplanetesimals is T 3 w 
9 • 10 3 h. At fc r =133, on the other hand, T s « 0.5 h. For 
greater values of a, T s is even smaller. Such small periods 
of axial rotation cannot exist, especially as bodies obtained 
by the contraction of rotating rarefied preplanetesimals are 
loosely bound, and can lose material easier than solid bodies. 

The value of v s = 2nrf/T s (the velocity of a particle 
on the surface of a rotating spherical object of radius 77 at 
the equator) is equal to v c f = (Gmf/rj) 1 ^ 2 (the minimum 
velocity of a particle that can leave the surface) at T sc = 
T s = (Stt/ P G) 1/2 . Thus, v s is equal to the equatorial escape 
velocity v c f at T se = T s — (3Tr/2pG) 1 ^ 2 . If k p is the density 



Considering that T sc = 2-wJ s /K s e, the critical radius at 
which material begins to leave a contracting preplanetesimal 
is equal to 



0.16x 2 m^G 



k cr k.Q kjjTu 
X 2 ' 



(12) 



in g cm' 3 , then T sc ~ 3.3/k 1 / 2 h and T se w 2.33/fcp 7 ^ h. 



.1/2 



The coefficient k cr in equation (12) is approximately 0.82 at 
mi = m-2 = raj /2, and k cr ~ 2.08(m2/mi) 2 at m2 mi w 
m/. Therefore, at fee = fc^ = x = 1 and mi = 7712, the ra- 
dius of the disc formed around the primary could be as much 
as 0.82rj^. In most cases, however, the disc would be rela- 
tively small since the masses of colliding preplanetesimals 
are usually different and the preplanetesimals are smaller 
than their Hill spheres. 

If preplanetesimals are smaller than their Hill spheres, 
it may be better to consider the encounter of the preplan- 
etesimals up to the radius of the Hill sphere corresponding 
to the mass rrif = mi + m-2, but not up to the sum of the 
radii of their Hill spheres, as it was in formulae (2), (5) and 
(10). For two identical preplanetesimals, this consideration 
decreases K a e and increases T a by a factor of 4 2//3 « 2.5. 



3.4 The spin imparted to a preplanetesimal by 
the accretion of much smaller objects 

A preplanetesimal can change its spin and mass due to the 
capture of much smaller objects (e.g. dust and boulders) ini- 
tially moved outside of its Hill sphere. In this subsection, we 
consider only accreted objects initially moving in circular he- 
liocentric orbits. At Aa = O ri2, the difference in the veloci- 
ties of two objects is about 0.5w c Ori2/a, where ri2 = ri+r 2 - 
The relative distance covered with this velocity during 
one revolution around the Sun is (0.5O« c ri2/a)(27ra/?; c ) = 
7rOri2. Thus, about (27ra)/(7rOri) = 2a/(Ori) revolutions 
are needed for an object with radius n to sweep up all 
smaller objects with radii r 2 <C ri. For example, taking 
Pluto's Hill sphere and semimajor axis for n and a, we 
obtain 2o/ri « 1.5 • 10 4 . This value may be similar to the 
contraction times of trans-Neptunian preplanetesimals, or 
smaller by one or two orders of magnitude (see the discus- 
sion in Section 1.1). 

A cross-section of a sphere is a circle, and the length of a 
chord at distance 6 n from its centre is 2n(l — O 2 ) 1 / 2 . The 
relative velocity at this distance is proportional to Qri/a. 
Therefore, if small objects are uniformly distributed around 
the orbit of the first object, the number of the small objects 
captured by a preplanetesimal per unit of time at Aa = 6n 
is proportional to (1 — 2 ) ly ^ 2 0. Multiplying this value by 
fc e (as K s oc fee), we obtain K s oc 9(1 - 6 2 ) 1/2 (1 - 1.59 2 ). 
By integrating over O from to 0.8165 and from 0.8165 to 
1 for a circular cross-section, I find that the ratio of positive 
to negative angular momenta is 9.4. 

Let us consider a spherical preplanetesimal of mass m pp 
and radius r — r pp that grows due to collisions with smaller 
objects as described above. Typical tangential velocity of 
collisions is supposed to be v T — 0.6v c r/a, and AK is the 
difference between fractions of positive and negative angu- 
lar momenta acquired by the preplanetesimal. The ratio be- 
tween the radii of a preplanetesimal of density p and the 
planetesimal formed by its contraction is denoted by k r . 
Considering an integral over a radius r of a growing preplan- 
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etesimal and supposing dK s — rv T dm and dm = 4np r 2 dr, 
we can calculate K s . 

Finally, the rotation period of the planetesimal formed 
by contraction of the preplanetesimal is 

T s = 2nJ s /K 3 « 7xa 3/2 (GMs)~ 1/2 AK- 1 k- 2 . 

If r pp is equal to the radius of the Hill sphere, we have 

T s = 5.7 X Ms /6 (Gay 1/2 AK- 1 p- 2/3 kr 2 - 

If the preplanetesimal is much smaller than its Hill sphere 
and v T — a v par (where v par is the parabolic velocity at the 
surface of the preplanetesimal), then 

K s « 0.85a(r pp G) 1/2 m 3 pp 2 AK « 0.67aG 1/2 P a 1/a m% 3 AK 
and 

T a « lA5 X (aAK)- 1 (p G)- 1/2 k- 2 . 

Based on simulations of a large number of encounters 
between two objects of density p moving around the Sun, we 
calculated values of AK for different elements of heliocen- 
tric orbits. The calculations showed that the fraction K+ 
of positive angular momenta (note that AK = 2K + — 1) 
tends to decrease with increases in a, p, or eccentricity. For 
RPPs and small eccentricities e, positive angular momen- 
tum increments are more frequent than negative increments 
(K+ > 0.9 at e=0). For solid objects in orbits with a = 1 au 
and e-(mi/Ms) -1 ' 3 ~ 2 — 7, AK is positive and about a few 
hundredths. For circular orbits, a was typically calculated 
to be 0.6. 



data used in these calculations were taken from the site 
|http:/ /www .joh nstonsarchive.net /astro , which also provides 
many references for each system. Similar (but older) data 
can be found in Richardson & Walsh (2006). For calcula- 
tion of tk ps = (Kscm + K sp in)/K 3 oe P s and rKeq = (K 3crn + 
Ks P in)/K s o6eq, only the spin of the primary was considered. 
However, these ratios would be about the same if they in- 
cluded the spins of all components, because K sp in is several 
times smaller than K SC m even for equal masses. For unequal 
masses most of the spin momentum is due to the primary 
component. 

Table 1 also reports the values of 2L/d p and L/vh, 
where th is the Hill radius for the total mass m ps of the 
binary. Three velocities are presented in the last lines of Ta- 
ble 1: v Tpr oe is the tangential velocity v T of the encounter 
at fee = 0.6 between Hill spheres with masses equivalent to 
the present primary and secondary components; t> T e 3 06 is the 
same velocity for Hill spheres of equal masses (0.5m ps ); and 
v eBC - P r is the escape velocity at the edge of the Hill sphere 
of the primary. 

A relative velocity v re i equal to 1 m s _1 is obtained for 
two colliding objects orbiting at e ~ 0.0002 and a=40 au. 
The model described above uses perfectly circular heliocen- 
tric orbits, and may not be correct at e > 0.0002 for preplan- 
etesimals corresponding to solid objects of diameter d ^ 200 
km because v re i ~ v e3C - pr for e ~ 0.0002 and d ~ 200 km. 
This conclusion is valid for other values of a, since both v re i 
and proportional to a 1//2 . It is reasonable to as- 

sume that the eccentricities of rarefied preplanetesimals are 
very small because of their interaction with gas and dust. 



4 COMPARISON OF ANGULAR MOMENTA 
OF DISCOVERED BINARIES WITH MODEL 
ANGULAR MOMENTA 

4.1 Description of the data presented in Table 1 

For six discovered binaries, the angular momenta K 3crn of 
the systems were estimated and compared to the present 
model. The primary and secondary components of the bi- 
nary systems have diameters d p and d s and masses m p and 
m s , reported in the first rows of Table 1. Using the formulae 
presented in Sections 3.2 and 3.3, I also calculate the the- 
oretical angular momentum K s o6 P s of two preplanetesimals 
with masses m p and m B encountered from circular helio- 
centric orbits at the mean value of fee (equal to 0.6), and 
the angular momentum K s o6eq of two preplanetesimals with 
masses equal to 0.5(m p + m s ). Thus, each system is com- 
pared to a theoretical scenario where two rarefied preplan- 
etesimals on different orbits merged then contracted into a 
pair of solid bodies. All three momenta are relative to the 
centre of mass of the system. 

While K s o6 P s and K a oeeq are calculated for RPP diam- 
eters equal to the Hill radius, binary systems can also be 
formed by the collision of much smaller preplanetesimals, as 
discussed in Section 3.3. Naturally, the heliocentric orbits of 
two preplanetesimals can have other separations, so fee can 
be smaller or larger than 0.6. 

The spin momentum of the primary is K ap i„ = 
0.27T x m P d p /Ts P , where T sp is the period of axial ro- 
tation of the primary. The separation distance between 
the primary and the secondary is denoted by L. The 



4.2 Discussion of the data presented in Table 1 

For the trans-Neptunian binaries presented in Table 1, nei- 
ther TK P s = {Kscm + Kspin) / Ks0f >p s nor VKeq = [K scm + 
K spin) I K sOQeq \T Kps ^ rtceq) exceed 0.63. This highest value 
was obtained for 2001 QW322, which consists of two approx- 
imately equal masses separated by a large distance. For the 
other trans-Neptunian binaries, L/th and, therefore, rpceq 
are smaller (often much smaller) than those for 2001 QW322- 
The small values of rKpa and rKeq can be explained if the 
discovered binaries arose from preplanetesimals that were al- 
ready partly compressed (i.e. smaller than their Hill spheres) 
and/or denser towards the centre at the moment of collision. 
Petit et al. (2008) noted that most models of binary forma- 
tion cannot easily explain the properties of 2001 QW322- For 
this binary, the equation K a & = K scm is fulfilled if fee ~ 0.4 
and v T ~ 0.16 m s^ 1 . Therefore, formation of this binary 
system can be explained as the merger of two RPPs in cir- 
cular heliocentric orbits. 

Another property of 2001 QW322 is that its components 
have a retrograde mutual orbit. As shown in Section 3.2, the 
angular momentum is negative for encounters of RPPs at 
0.82 < Q < 1. This condition can also result in the formation 
of widely separated binaries. 

Main-belt asteroid (87) Sylvia (384x262x232 km) has 
two satellites, Romulus (Li~1356 km, di«18 km) and Re- 
mus (I/2«706 km, d2~7 km). For this system, rueq ~ 
0.05 <H VKpa ~ 130. If the Sylvia system originated in the 
collision of two rarefied preplanetesimals, then their masses 
did not differ by more than an order of magnitude (as 

K sOGeq I K spin ~ 20). 
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Table 1. Angular momenta of several binaries. 



binary 


Pluto 


(90842) Orcus 


2000 CFi 05 


2001 QW322 


(87) Sylvia 


(90) Antiopc 


a, au 




oy.o 


A1 8 


A1 QA 




^ 1 

0. 1DO 


dp , km 




you 


1 7n 
1 / u 


1UO . 


ZOO 


oo 


d s , km 


1212 


260 


120 


108? 


18 


84 


m P 1 


1.3 x 10 22 


7.5 x 10 20 


2.6 x 10 18 ? 


6.5 x 10 17 ? 


1.478 x 10 19 


4.5 x 10 17 


777 „ Ircr 


1.52 x 10 21 


1.4 x 10 19 


9 x 10 17 ? 


6.5 x 10 17 ? 


3 x 10 15 


3.8 x 10 17 






tVkV ri — 1 cr f*m ^ 

1U1 iJ _1 • \J cL l_.ll! 






fr\r / 1 — 1 cf r*TYi 3 

IVJl /_/ _L ^.111 




km 


19,750 


8700 


23,000 


120,000 


1356 


171 


L/r H 


0.0025 


0.0029 


0.04 


0.3 


0.019 


0.007 


2L/dp 


16.9 


18.3 


271 


2200 


9.5 


3.9 


T S p, h 


153.3 


10 






5.18 


16.5 


-ffscm, kg km 2 s" 1 


6 x 10 24 


9 x 10 21 


5 x 10 19 


3.3 x 10 19 


10 17 


6.4 x 10 17 


K spin , kg km 2 s _1 


10 23 


10 22 


1.6 x 10 18 


2 x 10 17 


4 x 10 19 


3.6 x 10 16 








at T s =8 h 


at T s =8 h 






^ s06ps, kg km 2 s" 1 


8.4 x 10 25 


9 x 10 22 


1.5 x 10 20 


5.2 x 10 19 


3 x 10 17 


6.6 x 10 18 


K „06eq, kg km 2 s _1 


2.8 x 10 26 


2 x 10 24 


2.7 x 10 20 


5.2 x 10 19 


8 x 10 20 


6.6 x 10 18 


(^scm ~T" K spin) / K sQGps 


0.07 


0.2 


0.3 


0.63 


130 


0.1 


(Kscm ~\~ K spin} / K sOGeq 


0.02 


0.01 


0.2 


0.63 


0.05 


0.1 


VreqOI}, m S _1 


6.1 


2.2 


0.36 


0.26 


2.0 


0.82 


tVpr06, m S _1 


5.5 


1.8 


0.3 


0.26 


1.3 


0.82 


Vesc — pr , m S 


15.0 


5.8 


0.8 


0.53 


5.3 


1.7 



Main-belt asteroid (90) Antiope belongs to a syn- 
chronous double-asteroid system. Its rotational and orbital 
periods are 16.5 h, and the mutual orbit of the components is 
almost circular. Before tidal forces induced these orbits, the 
components may have been separated by a distance even 
smaller than the present 171 km. Since the present ratio 
2L/d v is small, this binary likely formed at the stage of solid 
(or almost solid) bodies. For this binary, tk P s ~ TKeq ~ 0.1. 
Since K SC m/ 1 K s e oc (L/a) 1 ^ 2 (see Section 3.3), the model of 
rarefied Hill sphere preplanetesimals in almost circular helio- 
centric orbits is capable of producing a binary system with 
masses of components similar to those for Antiope. However, 
the separation L obtained is greater by two orders of mag- 
nitude than that of the Antiope binary and is comparable 
to the radius of the Hill sphere. 

4.3 Formation of Pluto's axial rotation and the 
inclined mutual orbits of binary systems 

Pluto has three satellites, but the contribution of two satel- 
lites (other than Charon) to the total angular momentum of 
the system is small; for the whole system K SC m/ 'K sp in ~ 60. 
The axial tilt of Pluto with respect to its orbit is 119.6° . Such 
a reverse rotation is possible in the collision of two (solid or 
rarefied) objects, but at the stage of rarefied preplanetesi- 
mals, it is better to use | fee | ~ 0.2 — 0.3 (see Section 3.2). 
In the proposed model, it is not possible to obtain Pluto's or 
any reverse rotation from a large number of collisions with 
smaller objects (as in Section 3.4). 

For a collision of two identical Hill sphere preplanetes- 
imals, the value of K s & obtained for | fee | = 0.3 exceeds 
Kscm for the Pluto system by a factor of 23. If the masses of 
the parent Hill sphere RPPs are equal to the masses of Pluto 
and Charon, then K a e/Ks C m ~ 7 at | fee | = 0.3. Further- 
more, the angular momentum of colliding preplanetesimals 
in eccentric orbits can be greater than that obtained for cir- 
cular orbits. The above estimates imply that the radii of the 
preplanetesimals that gave rise to the Pluto system were 



smaller by at least an order of magnitude than their Hill 
radii. 

Most of the angular momentum of the Pluto system 
could have resulted from a single collision of preplanetesi- 
mals moving in different planes. To explain Plutos axial tilt 
(about 120°) and the inclined mutual orbit of 2001 QW322 
(124° from the ecliptic) in these terms, I first note that the 
thickness of the disc (i.e., the distances between the centres 
of mass of the colliding preplanetesimals and the equato- 
rial plane of the disc) was comparable to or greater than 
the radii of the RPPs themselves. The inclined mutual or- 
bits of many trans-Neptunian binaries support the idea that 
the momenta of such binaries mainly resulted from single 
collisions of preplanetesimals rather than the accretion of 
smaller objects. (In the latter case, the primordial inclina- 
tions of mutual orbits to the ecliptic would be relatively 
small.) 

In the models of binary formation considered here, the 
spin vector of the primary preplanetesimal is almost perpen- 
dicular to the plane in which satellites of the primary move 
(as is the case with Pluto). However, the direction of this 
vector could vary over time due to collisions of the primary 
with solid bodies. 



5 DISCUSSION 

5.1 Radii and mergers of rarefied 
preplanetesimals 

The radii of preplanetesimals need to be comparable to their 
Hill radii only if the resulting binary is to have a separation 
L comparable to its own Hill radius th- Excepting 2001 
Q W322 , all the binaries presented in Table 1 and many other 
discovered binaries have ratios L/rn below 0.04. Among the 
binaries considered by Richardson & Walsh (2006) in their 
fig. 1, this ratio was less than 0.3. RPPs much smaller (by at 
least by an order of magnitude) than th suffice to explain the 
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formation of binaries with L/rn < 0.04. Generally speaking, 
the required radii of the rarefied preplanetesimals can be of 
the same order as the present distance L between the binary 
components. However, the secondary component could form 
much closer to the primary than the outer edge of the disc, 
and L could increase due to tidal forces. This estimate of the 
preplanetesimal size is therefore not an accurate predictor. 

The values of 2L/d p presented in Table 1 vary from 17 
to 2200 for trans-Neptunian binaries. The trans-Neptunian 
binaries considered by Richardson & Walsh (2006) fall inside 
the same range, but asteroid binaries may have smaller val- 
ues. Cuzzi et al. (2008) considered spherical, rarefied clumps 
of diameter I = (1 — 5) x 10 4 km and mass equivalent to 
a body of unit density and radius 10-100 km, orbiting at 
a=2.5 au. The diameters I are greater than the separations 
of known asteroid binaries. Rarefied pre-asteroids may there- 
fore have solidified before colliding with other pre-asteroids. 

The densities of rarefied preplanetesimals can be very 
low, but their relative velocities v re i during an encounter 
are also very small. In particular, the relative speed of the 
collision is typically smaller than the escape velocity v e3C - pr 
at the Hill radius of the primary. In addition to v esc - pr , Ta- 
ble 1 presents the tangential component v Tpr oa of v re i- In 
some models of the evolution of Saturn's rings (R. Perrine, 
private communication, 2009) colliding objects form a new 
object if their impact speed is less than the mutual escape 
speed by a certain factor. If the greater encountered object 
is much smaller than its Hill sphere, and if both heliocentric 
orbits are almost circular before the collision, then the ve- 
locity of the collision v co i w {v 2 e i + n^) 1 ' 2 does not differ 
much from the parabolic velocity v par at the surface of the 
primary RPP (radius r pc ). Indeed, v par is proportional to 
r P c^ 2 . Therefore, encounters could result in a merger (fol- 
lowed by formation of a binary) at any r pc < th- 

For some pairs of objects, the mutual orbit can be more 
complicated than a parabola or ellipse. In a coplanar model, 
decreasing r pc by a factor of 10 can reduce the fraction 
of encounters resulting in a collision by a factor less than 
10. Greenzweig and Lissauer (1990) published plots show- 
ing the dependence of the closest approach r min between 
objects in initially circular heliocentric orbits separated by 
an initial distance Aa. Based on these figures, the maxi- 
mum value of Aa is 0A7th for r min < Q.lrH and Q.lQrn 
for r min < 0.01r£f. In Greenzweig and Lissauer's study, the 
integration began at a distance greater than th and consid- 
ered close encounters taking place at Aa < 3.5vh- Note that 
in the above example, while the densities of uniform spheri- 
cal preplanetesimals ranging from vh to O.Olr// in size differ 
by six orders of magnitude, the probability of collision dif- 
fers only by a factor of ~ 20 provided all motions occur in 
the same plane. The difference in the probabilities is much 
greater in a non-planar model. 

Johansen et al. (2007) determined that the mean free 
path of a boulder inside a cluster - preplanetesimal is shorter 
than the size of the cluster. This result supports the pic- 
ture of mergers between rarefied preplanetesimals. Lyra et 
al. (2008) noted that the velocity dispersion of rarefied pre- 
planetesimals remains below 1 m s _1 in most simulations, 
so destructive collisions between boulders are avoided. 

To illustrate the probability of a merger between two 
rarefied preplanetesimals, consider the following simple 
model. There is a spherical preplanetesimal of diameter D s 



and mass M, consisting of N identical boulders of diameter 
d. After contraction, a solid, spherical planetesimal of diam- 
eter D is formed. The densities p of the boulders and the 
solid planetesimal are the same. A second preplanetesimal 
then passes through this cluster. The ratio of the length of 
its path inside the first preplanetesimal to D s is denoted as 
k s . Because the relative motion of the centres of mass of the 
preplanetesimals can be complicated, k s can exceed 1. 

The boulders belonging to both preplanetesimals are 
considered to have identical diameters d, and the volume 
swept by one boulder is -k d 2 k s D s . The ratio of this volume to 
the volume ^D 3 of the Hill sphere, divided by N = (D/d) 3 , 
gives the number of collisions N co i — 6k s D 3 /(D^d). Let us 
take D s = 2k H a(M/3M s ) 1/3 , p = k p g cm~ 3 , and a = k a 
au. In this case, we obtain 

D/d w N co i x 3 • 10 3 ■ klk 2 / 3 ^ 1 ^. (13) 

For D = 1000 km, d = 0.3 m, and k p = k 3 = k H = N co i = 1, 
equation (13) is fulfilled at a — 33 au. Using equation (13), 
we find that N co i oc D ■ k^ 2 . This verifies the intution that 
for small values of the ratio kn between a preplanetesimal's 
radius and its Hill radius, N co i can be relatively large. N co i 
also decreases with D, but at D = 50 km and kn = 0.2, 
Ncoi is almost the same as for D — 1000 km and kn = 1- 
This means that for most binaries in the Solar system with 
D > 50 km, their parent preplanetesimals should have had 
N co i ^ 1 provided that d $C 1 m. Boulders (or dust particles) 
are more likely to be captured if they have smaller diame- 
ters d. The probability of capture also increases for boulders 
(particles) closer to the centres of the preplanetesimals, if 
the density of the preplanetesimal is higher at a smaller dis- 
tance from the centre. The relative velocities of the boulders 
in this model are usually smaller than the escape velocities 
of the preplanetesimals, so some collided boulders could re- 
main inside the Hill sphere. 

At greater eccentricities, the mean collision velocity be- 
tween preplanetesimals and the minimum distance of closest 
approach between material points corresponding to encoun- 
tering preplanetesimals are greater, and the particles remain 
inside the Hill sphere for less time. Therefore, the probabil- 
ity of a merger of preplanetesimals is smaller. However, the 
typical angular momentum of preplanetesimals encountered 
up to the Hill sphere is greater. 

The fraction of trans-Neptunian binaries that formed 
during the preplanetesimal stage depends on the initial dis- 
tribution of sizes, densities, and contraction times among 
preplanetesimals. These issues must wait on further studies 
of the formation and evolution of preplanetesimals. 



5.2 Binaries that originated at different distances 
from the Sun 

Given a primary of mass m p and a much smaller sec- 
ondary, both in circular heliocentric orbits, one can obtain 
v T /v e3 c- P r ~ 0.66fee(-^ £ ) 1 ^ 3 /a (see designations and for- 
mulae in Sections 3.2 and 4.1). This ratio decreases with 
increasing a and m p . Therefore, preplanetesimals are more 
likely to merge when the primary is more massive and lo- 
cated farther from the Sun. The total mass of all preplan- 
etesimals and the ratio of the time needed for contraction of 
preplanetesimals to the period of their rotation around the 
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Sun could be greater for the trans-Neptunian region than 
for the initial asteroid belt. (Several authors have arrived at 
the similar conclusion on the ratio for dust condensations, 
e.g. Safronov 1969.) 

The above factors could explain why a larger fraction 
of binaries are found at greater distances from the Sun, and 
why the typical mass ratio (secondary to primary) is greater 
for TNOs than for asteroids. The binary fractions in the mi- 
nor planet population are about 2 per cent for main-belt 
asteroids, 22 per cent for cold classical TNOs, and 5.5 per 
cent for all other TNOs (Noll 2006). Note that TNOs mov- 
ing in eccentric orbits (the third category) are thought to 
have formed near the giant planets, closer to the Sun than 
classical TNOs (e.g. Ipatov 1987). 

5.3 Asteroid binaries 

In our opinion, some asteroid binaries with large primaries 
could have formed at the preplanetesimal stage. The eccen- 
tricities of satellite orbits around large asteroids (including 
Sylvia and Antiope) are usually (but not always) relatively 
small. Under the proposed model of rarefied preplanetesimal 
mergers, small eccentricities could indicate that the satel- 
lites formed from a disc around the primary. The plane of 
the disc will be close to that defined by the line connect- 
ing two preplanetesimals and their relative velocity vector. 
If the initial preplanetesimal orbits have large eccentricities 
e and/or inclinations i, the relative velocity depends on the 
angle between the planes of heliocentric orbits. For small e 
and i, the plane of relative motion inside the Hill sphere can 
differ significantly from the orbital planes. 

If the thickness of the disc of preplanetesimals is of the 
same order as the diameters of the largest Hill spheres, then 
the proposed model of binary formation permits any incli- 
nation i m between the mutual orbit of a binary and the 
ecliptic. Some asteroid binaries have highly inclined mutual 
orbits (e.g. i m —64° for Antiope and i m =93° for Kalliope). 
For (87) Sylvia, (107) Camilla, and (283) Emma, on the 
other hand, i m does not exceed 3°. High- velocity collisions 
between solid asteroids may more often yield i m of the or- 
der of typical inclinations in the asteroid belt, as most of the 
debris will have neither large nor very small inclinations. 

The total angular momentum of two identical Hill 
sphere preplanetesimals in circular heliocentric orbits often 
exceeds (In the case of (87) Sylvia, by an order of magni- 
tude) the critical value at which particles leave the surface 
of the solid primary. 

Solid bodies can also attain the critical angular mo- 
mentum (corresponding to T s ^ SS/k^ 2 h; see the end of 
Section 3.3) by collisions, but eccentric (and/or inclined) 
heliocentric orbits are required. Using equation (6) with 
mi = 77i2, v T /x = 3.5 km s -1 (e.g. at a typical velocity 
v co i of collisions in the asteroid belt equal to 5 km s _1 , 
v T « 0.7v co i, and \ — 1), and 77 ~ 6600 km, I obtain 
T s ~ 3.3 h. Therefore, the critical angular momentum can 
also result from the collision of two identical solid asteroids 
of any radius (<6600 km). 

For a given ratio mi I mi and eccentricities of heliocen- 
tric orbits of solid bodies, the probability of attaining the 
critical momentum during a collision is greater for smaller 
values of mi and a (because T s oc rf/v T , and at smaller 
a both orbital velocity and v T are larger). For v r /x = 3.5 



km s _1 , T s = 3.3 h, and n = 100 km; formula (7) then 
yields mz/mi « 0.006. 

Solid bodies can be disrupted during a collision, whether 
or not they reach the critical angular momentum. Most of 
the material ejected in this fashion would leave the Hill 
sphere, especially when the primary has a relatively small 
mass. The spin periods of asteroid primaries with d > 50 
km are often several times greater than 3 h. It may be that 
their original rotation was faster. 

Some present-day asteroids (especially those with d < 
10 km) may be debris from larger solid bodies. As discussed 
in Section 1.2, binary asteroids in the near-Earth object pop- 
ulation formed at the stage of solid bodies. The spin and 
shape of these bodies could have changed during the evolu- 
tion of the Solar system due to collisions with other small 
solid bodies. 

5.4 The colours and total mass of TNOs 

In the proposed model of binary formation, two preplanetes- 
imals originate at almost the same distance from the Sun. 
This point agrees with the correlation between the colours 
of primaries and secondaries obtained by Benecchi et al. 
(2009) for trans-Neptunian binaries. In addition, the ma- 
terial within the preplanetesimals could have been mixed 
before the binary components formed. 

The formation of classical TNOs from rarefied preplan- 
etesimals could have taken place for a small total mass of 
preplanetesimals in the trans-Neptunian region, even given 
the present total mass of TNOs. Models of TNO formation 
from solid planetesimals (e.g. Stern 1995; Davis & Farinella 
1997; Kenyon & Luu 1998, 1999) require a more massive pri- 
mordial belt and small (~ 0.001) eccentricities during the 
process of accumulation. However, the gravitational inter- 
actions between planetesimals during this stage could have 
increased the eccentricities to values far greater than those 
mentioned above (e.g. Ipatov 1988, 2007). This increase tes- 
tifies in favour of formation of TNOs from rarefied preplan- 
etesimals. 



6 CONCLUSIONS 

This analysis has shown that some trans-Neptunian objects 
could have acquired their axial momenta and/or satellites 
during a primordial stage as rarefied preplanetesimals. In 
this scenario, most rarefied pre-asteroids could have solidi- 
fied before colliding with other pre-asteroids. 

Models of binary formation due to gravitational inter- 
actions or collisions between objects within the Hill radius, 
which have been studied by several authors for solid objects, 
could be more effective for rarefied preplanetesimals. For ex- 
ample, due to their almost circular heliocentric orbits, two 
nearby preplanetesimals remain inside the Hill radius longer. 
The centres of mass of two rarefied preplanetesimals may 
also be able to approach each other more closely than solid 
planetesimals because the planetesimals could have greater 
eccentricities of heliocentric orbits. 

During a collision, rarefied preplanetesimals can have 
higher densities closer to their centres. In this case, there 
could be two centres of contraction inside the rotating pre- 
planetesimal formed as a result of the collision of two rar- 
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efied preplanetesimals. The result would be a binary with 
two roughly equal masses, which could be separated by any 
distance up to the Hill radius. The eccentricity and inclina- 
tion of the secondary component's orbit around the primary 
component can have any value. The observed separation dis- 
tance can be of the order of the radius of a greater encoun- 
tered preplanetesimal. 

Some binaries could form because the angular momen- 
tum of a binary that was obtained at the stage of rarefied 
preplanetesimals was greater than the angular momentum 
that can exist for solid bodies. Material that left a contracted 
preplanetesimal formed as a result of a collision of two pre- 
planetesimals could form a disc around the primary. One 
or more satellites could grow in the disc within the Hill ra- 
dius, typically much less. These satellites move mainly in 
low-eccentricity orbits. 

The two scenarios described above can take place at 
the same time. It is thereby possible that, besides massive 
primary and secondary components, smaller satellites can 
form around the primary and/or the secondary. 

Among discovered trans-Neptunian binaries, the an- 
gular momentum is usually considerably smaller than the 
model's prediction for two identical rarefied preplanetesi- 
mals having the same total mass and encountering up to 
the Hill radius from circular heliocentric orbits. The predic- 
tions for two preplanetesimals with unequal masses equiva- 
lent to those observed in the trans-Neptunian binaries are 
also too large. Furthermore, the observed separations be- 
tween components are usually much smaller than the Hill 
radii. These facts support the hypothesis that most preplan- 
etesimals were already partly compressed at the moment of 
collision, i.e., smaller than their Hill radii and/or centrally 
concentrated. The contraction of preplanetesimals could be 
slower farther from the Sun, which would explain the greater 
fraction of binaries formed at greater distances from the Sun. 

The results of this research, which has focussed on the 
angular momentum of colliding preplanetesimals, can also 
be used to analyse the formation of axial rotation of rarefied 
and solid bodies. 
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